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Multichannel Quantum Defect Theory (MQDT) is shown to be capable of producing quantita- 
tively accurate results for low-energy atom-molecule scattering calculations. With a suitable choice 
of reference potential and short-range matching distance, it is possible to define a matrix that en- 
capsulates the short-range collision dynamics and is only weakly dependent on energy and magnetic 
field. Once this has been produced, calculations at additional energies and fields can be performed 
at a computational cost that is proportional to the number of channels TV and not to A*'^. MQDT 
thus provides a promising method for carrying out low-energy molecular scattering calculations on 
systems where full exploration of the energy- and field-dependence is currently impractical. 
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Note for copy editor: We have been very careful to 
make correct use of Roman and italic subscripts and 
superscripts, with Roman for abbreviations and italic 
for mathematical indices. Please do not change all our 
subscripts and superscripts to italic. Also note that 
subscript lower-case o is an abbreviation for "open" 
and should not be changed to zero. 



I. INTRODUCTION 

The creation of the first dilute atomic Bose-Einstein 
condensates (BECs) in 1995 0, Q led to enormous ad- 
vances in ultracold atomic physics. There is now great 
interest in producing samples of cold molecules, at tem- 
peratures below 1 K p-Q, and ultracold molecules, at 
temperatures below 1 mK There are many poten- 

tial applications of ultracold molecular samples, amongst 
which are high-precision measurements [TolfTTI , quantum 
computation ^1^] and ultracold chemistry |13| . 

Understanding atomic and molecular interactions and 
collisions is essential to the study of cold and ultra- 
cold molecules. For example, methods such as buffer- 
gas cooling and Stark deceleration [l^ can produce 
cold molecules with temperatures between 10 mK and 
1 K. However, a second-stage cooling method is needed 
to bring the molecules into the ultracold regime. Sym- 
pathetic cooling, in which the molecules are allowed to 
thcrmalize with a gas of ultracold atoms, is a promising 
second-stage cooling method [l^. However, while elas- 
tic collisions allow thermalization, inelastic collisions can 
cause trap loss 17[, and for many systems the inelastic 



collisions are predicted to be too large for sympathetic 
cooling to succeed [T8l - [20| . Scattering calculations are 
essential in order to identify systems for which sympa- 
thetic cooling has a good prospect of success. Once in the 
ultracold regime, the extent to which atomic and molec- 
ular interactions can be controlled again depends on a 
detailed understanding of their coUisional properties. 



Quantum molecular scattering calculations are usu- 
ally carried out using the coupled-channel method: the 
Schrodinger equation for scattering is converted into a set 
of coupled differential equations, which are then propa- 
gated across a range of values of the intermolecular dis- 
tance r. The size of the problem is determined by the 
number of channels N (the number of coupled equations). 
The usual algorithms take a time proportional to N^, 
since each step of the propagation requires an 0{N^) 
matrix operation. 

Cold molecule scattering presents problems with a 
large number of channels for two reasons: 

1. At very low energies, small splittings between 
molecular energy levels become important. This 
makes it necessary to include fine details of molec- 
ular energy level patterns, such as tunneling and 
nuclear hyperfine splitting. The extra degrees of 
freedom require additional basis functions; in par- 
ticular, including nuclear spins can multiply the 
number of equations by a substantial factor (some- 
times 100 or more). 

2. Collisions in the presence of electric and magnetic 
fields are very important. In an applied field, the 
total angular momentum J is no longer a good 
quantum number. Because of this, the large sets 
of coupled equations can no longer be factorized 
neatly into smaller blocks for each J as is possible 
in field-free scattering. 

In addition, in cold molecule applications it is often neces- 
sary to repeat scattering calculations on a fine grid of en- 
ergies and/or applied electric and magnetic fields, which 
adds greatly to the computational expense. 

Multichannel Quantum Defect Theory (MQDT) of- 
fers an alternative to full coupled-channel calculations. 
It was originally developed to provide a uniform treat- 
ment of bound and scattering states for problems involv- 
ing the interaction an electron with an ion core with 
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Coulomb forces at long range [2l|, |22|, but was subse- 
quently generalized to handle a range of other long-range 
potentials [22l - [27j . It has been successfully applied to 
scattering problems as diverse as negative ion photode- 
tachment 1281. near-threshold predissociation of diatomic 
molecules |25l . [29l | and predissociation of atom-diatom 
Van der Waals complexes [13, HH . More recently it has 
been applied to ultracold collisions between pair s of neu- 
tral atoms [33 - l36| , between atoms and ions |37l . [ssj , and 
between highly reactive molecules [39l - l4l| . 



MQDT can be viewed in two different ways. The first 
tries to capture the important physics of collisions within 
a few analytic quantum defect parameters. The other 
views it as a method for solving the coupled equations 
of scattering theory which offers substantial insights and 
advantages in efficiency. The common feature of the two 
approaches is to take advantage of the enormous differ- 
ence in energy and length scales associated with sepa- 
rated collision partners and short-range potentials. 



II. THEORY 
A. Coupled-channel method 

Cold atomic and molecular collisions and near- 
threshold bound states are conveniently described by a 
set of coupled equations. The Hamiltonian for an inter- 
acting pair of atoms or molecules is of the form 



2^ 



V^ + Hint{T) + Vir,T), 



(1) 



where ^ is the reduced mass, is the Laplacian for the 
intermolecular coordinates, and r denotes all coordinates 
except the interparticle distance r. Hi^t^r) represents the 
internal Hamiltonians of the two particles and V{r,T) 
is the interaction potential. The total wavefunction is 
expanded 



5'(r,T) 



(2) 



When MQDT is viewed as a numerical method for solv- 
ing the coupled differential equatioiis, the goal is to ob- 
tain a matrix Y{E,B) [H, ^ [H, HI] that completely 
describes the short-range dynamics and is insensitive to 
collision energy E and magnetic field B. This matrix can 
be obtained once and then used for calculations over a 
wide range of energies and fields, or obtained by interpo- 
lation from a few points. MQDT achieves this by defin- 
ing Y{E, B) at relatively short range, as described below. 
The threshold behavior is accounted for from properties 
of single channels. Once the matrix Y{E, B) has been 
obtained, the time required for calculations at additional 
energies and fields is only proportional to iV, not . 



Understanding threshold atomic physics in quantum 
defect terms is weU developed [H^sllil El- Thresh- 
old bound-state and scattering properties are determined 
mainly by the long-range potential, which can often be 
approximated as — C„/r". For the case of the Van der 
Waals interaction, —Cq/t^, the linearly independent pair 
of solutions for a single potential is known [44!| . An ana- 
lytic approach to MQDT using these solutions has been 
developed jl^ an d g ives much insight into ultracold 
atom-atom collisions [43 . 



This paper investigates the use of MQDT as a numer- 
ical method to study cold atom-molecule collisions. The 
structure of the paper is as follows. In Section II we give 
an overview of the theory of MQDT, sufficient to define 
notation. In Section III we apply MQDT to the proto- 
type system Mg-|-NH, and compare it with full coupled- 
channel calculations in order to establish what is required 
for it to give accurate results. In Section IV we present 
our conclusions and suggest directions for future work. 



where the N functions </5i(r) form a basis set for the 
motion in all coordinates, r, except the intermolecular 
distance, and V'i(^) is the radial wavefunction in chan- 
nel i. Substituting this expansion into the total time- 
independent Schrodinger equation and projecting onto 
the basis function (pj(r) yields the usual coupled equa- 
tions of scattering theory. 



2^1 dr^ 



E 



where E is the energy. The coupling matrix W has ele- 
ments 



Hi,,t{T)+V{r,T) 



+ 



2/ir^ 



ipi{T)dT, (4) 



where Li is the partial- wave quantum number for channel 
i. Equation (jSj can conveniently be written in matrix 
form, 



2^ dr2" 



[W{r) - EI]if){r), 



(5) 



where ip{r) is a column vector made up of the solutions 
TpiiT) and I is the identity matrix. 

For both bound states and collision calculations, the 
wavefunction must be regular at the origin. When 
V{r) ^ as r — ^ 0, the short-range boundary condition 
is 



as r — >• 0. 



(6) 



At any energy, there are N linearly independent solution 
vectors ip{r) that satisfy these boundary conditions, and 
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it is convenient to combine them to form the N x N 
wavefunction matrix 'S'(r). 

The coupled-channel approach propagates either the 
wavefunction matrix ^{r) and its derivative '^''(r), or the 
log-derivative matrix L{r) — outwards from 

r — (or a point in the deeply classically forbidden region 
at short range) [i^, H^. In scattering calculations, the 
propagation is continued to a point Tmax at large r. The 
wavefunction or log-derivative matrix is then transformed 
into a representation where W is asymptotically diagonal 
I50ll, such that 



h^Li{Li 



1) 



Sij + 1 



(7) 



where n is the power of the leading term in the poten- 
tial expansion and is the threshold of channel i. 
Each channel is either asymptotically open, E > E°° , or 
asymptotically closed, E < E°° . The scattering bound- 
ary conditions are 



* = J{r)+N{r)K. 



(8) 



The matrices J and N are diagonal matrices containing 
Riccati-Bessel functions for open channels and modified 
spherical Bessel functions for closed channels [49|. In a 
problem containing N channels. No of which are open, 
the scattering S matrix is related to the open-open sub- 
matrix of K by 



S = (1 + iKooy\l - ^Ko 



(9) 



In full coupled-channel calculations, the matrices K and 
S are rapidly changing functions of both energy and field, 
particularly near scattering resonances, so that the entire 
propagation to long range must be repeated for each set 
of conditions required. 



B. Multichannel Quantum Defect Theory 

MQDT also begins by propagating the wavefunction or 
log-derivative matrix outwards from short range. How- 
ever, instead of continuing to rmax, matching takes place 
at a point rmatch, at relatively short range. The match- 
ing in MQDT treats the open and weakly closed channels 
on an equal footing; weakly closed channels are usually 
defined as those that are locally open, E > Wu^r), at 
some value of r, so are capable of supporting scattering 
resonances. Matching at short range produces a matrix 
Y{E, B) that is relatively insensitive to energy and ap- 
plied field, as described below. Y also varies smoothly 
across thresholds, unlike 5' and K . Provided the chan- 
nels are uncoupled outside Tmatch, it is then possible to 
obtain the scattering S matrix from Y using the proper- 
ties of individual uncoupled channels. 

We consider a problem with No open channels and A^c 
weakly closed channels at some collision energy E and 
field B. For each such channel, i = 1,N, where iV = 



No + Nc, MQDT requires a reference potential, Ul^^{r), 
which asymptotically has similar behavior to Wii{r) in 
equation ([7]). This reference potential defines a linearly 
independent pair of reference functions fi{r) and gi{r), 



+ Kf{r) 



Mr) = 0, 



(10) 



and similarly for gi, where the local wave vector Ki{r) is 
.(r) = yp(i?-[/-f(r)). (11) 



The regular solution fi has the boundary condition 
/i — > as r ^ 0. /; and gi are normalized to have 
Wentzel-Kramers-Brillouin (WKB) form, with amplitude 
Ki{r)~'^/'^ ^ at some point in the classically allowed region 
[11]. The N X N matrix Y is defined by matching at 

^matchi 



^ = f{r)+g{T)Y, 
or in terms of the log-derivative matrix L, 
{Lf^f) = {Lg-g')Y, 



(12) 



(13) 



where / and g are diagonal matrices containing the func- 
tions fi and gi and the primes indicate radial derivatives. 

In order to relate Y to the physical scattering S ma- 
trix, the asymptotic forms of the reference functions /; 
and gi in each channel are required. To this end another 
pair of reference functions is defined for each channel. For 
open channels, these functions are asymptotically energy- 
normalized, 



s^{r) 
Q(r) ' 



fc, ^ sin fc,T 



k- ^ cos fc,T 



2 



(14) 
(15) 



where is the phase shift associated with reference po- 
tential i and ki is the asymptotic wave vector. 



h-xr^iE 



(16) 



These asymptotically normalized functions are related to 
fi and gi through the quantum defect parameters Ci and 
tan Xi , 

Mr) = C-'fdr); (17) 
Mr) = C,[5,(r)+tanA,/,(r)]. (18) 

Thus Ci relates the amplitudes of the energy-normalised 
reference functions to WKB-normalised ones, while 
tan Xi describes the modification in phase due to thresh- 
old effects. Far from threshold, Ci « 1 and tauAi « 0. 

For each weakly closed channel, an exponentially de- 
caying solution is defined. 



Mr) 



l„-\ki\r 



(19) 
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This is related to the solutions fi and gi by a normaliza- 
tion factor Afi , and an energy- dependent phase i/i , 



4)t{r) = Mi [cosuif.i{r) - sint'igi(^)] 



(20) 



The phase Vi is an integer multiple of tt at each energy 
that corresponds to a bound state of the reference poten- 
tial in channel i. 

The y matrix is converted into the S matrix of scat- 
tering theory using the quantum defect parameters Ci, 
tan Ai, iaiiVi and ^i. First, the effect of coupling to closed 
channels is accounted for, 



Y = Yoo - roc[tani/ + Y,,]-^Y,, 



(21) 



where tanf is a diagonal matrix of dimension x 
containing elements iaiWi. The No x A'o matrix Y in- 
corporates any resonance structure caused by coupling to 
closed channels through tani/. Unlike Y itself, Y can be 
a rapidly varying function of energy and field. Secondly, 
threshold effects from asymptotically open channels are 
incorporated. 



R = C 



— tan A 



C 



(22) 



where C and tan A are diagonal matrices of dimension 
No X No, containing elements Ci and tanA^. Finally, the 
S matrix is obtained from 



S = e'^ [l + iR\[l-iR\' 



(23) 



This may be compared to equation (|9]) for the full 
coupled-channel method. The inclusion of the diagonal 
matrix e*^ accounts for the phase difference between the 
reference functions fi and gi used by MQDT and the 
Riccati-Bessel functions used by the full coupled-channel 
method. 

The approach taken in the present paper is somewhat 
different from that in refs. [l^llJ]. There MQDT was ap- 
proached as an exact representation of the full coupled- 
channel solution. The matrix Y was evaluated at a dis- 
tance Tinatch large enough that it had become constant 
as a function of rmatch- When this is done, MQDT gives 
the same (exact) results for any choice of reference poten- 
tial U^°^{r), although constancy of Y may be achieved at 
different values of rmatch for different choices. In our ap- 
proach, r„iatch is chosen to ensure that Y is only weakly 
energy-dependent, and this may require matching in a 
region where Y is not yet independent of rmatch- With 
this approach, MQDT provides an approximate solution 
whose quality depends on the choice of reference poten- 
tials. 



C. Numerical evaluation of reference functions and 
quantum defect parameters 

1. Open channels 

For an open channel i, the reference function Si is ob- 
tained by propagating a regular solution of (jlO[) from a 



point inside rmatch to a point rmax at long range and im- 
posing the boundary condition (or its Bessel function 
equivalent). This establishes the normalization of Si and 
also gives the phase shift , which is then used to obtain 
the function Ci at rmax from the boundary condition (jisp . 
The reference function Ci is then propagated inwards to 
rmatch- The two remainiiig quantum defect parameters 
are obtained by applying [25| 



and 



cot A,; 



Ki{ji ~ Ui) 



liUi 



(24) 



(25) 



in the classically allowed region, where 7i = s[/ Si and 
Ui — c'J Ci. The primes indicate radial derivatives. Equa- 
tions pT|) and (IT8|) then give the reference functions fi 
and gi. 



2. Closed channels 

For a weakly closed channel i, the reference function 
/,; is again obtained by propagating a regular solution of 
(jlOP outwards from a point inside rmatch, but in this case 
fi is normalized in the classically allowed region such that 

fUKl+ih=K,. (26) 

In the closed-channel case, gi cannot be obtained directly 
from fi at a single point. Instead, the reference function 
(j)i is obtained by using (jl9p as a long-range boundary 
condition and propagating a solution of (jlOl) inwards to- 
wards r — 0. The normalization factor Mi of equation 
(I^Ul) is obtained by matching to 

M^ ^{4>^J<, + ct>f/K,). (27) 

in the classically allowed region. The quantum defect 
parameter tan Vi is then obtained from 



tanvi 



Ktjti - ji) 



(28) 



where ti — (l)[/(j)i. Finally, the function gi is obtained 
from fi and 0i using equation (|20p . 



D. Sources of Error 

There are a number of sources of errors in MQDT cal- 
culations using our approach: 

1. Interchannel couplings that occur outside rmatch, 
which are not taken into account by equations (|2ip 
to (ESI; 

2. Deviations between the reference potentials Ul^^{r) 
and Wii(r) outside rmatch; 

3. Differences between the actual Y matrix at a given 
energy and field and the Y matrix obtained by in- 
terpolation. 
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III. RESULTS AND DISCUSSION 

To explore the application of MQDT to cold molec- 
ular collisions, we consider the prototype system 
Mg-|-NH(^I]~). The potential energy surface for this sys- 
tem is moderately anisotropic [5l| and provides substan- 
tial coupling between channels. The system is topical 
because Wallis and Hutson [s^] have shown that sym- 
pathetic cooling of cold NH molecules by ultracold Mg 
atoms has a good prospect of success. 

The energy levels of NH in a magnetic field are most 
conveniently described using Hund's case (b), in which 
the molecular rotation n couples to the spin s to produce 
a total monomer angular momentum j. In zero field, 
each rotational level n is split into sublevels labeled by 
j. In a magnetic field, each sublevel splits further into 
2j + 1 levels labeled by m^, the projection of j onto the 
axis defined by the field. For the n = levels that are 
of most interest for cold molecule studies, there is only 
a single zero-field level with j — \ that splits into three 
components with rrij = and —1. 

The coupled equations are constructed in a partly cou- 
pled basis set \nsjmj)\LML) , where L is the end-over-end 
rotational angular momentum of the Mg atom and the 
NH molecule about one another and Ml is its projec- 
tion on the axis defined by the magnetic field. Hyperfine 
structure is neglected. The matrix elements of the total 
Hamiltonian in this basis are given in ref . (50| . The only 
good quantum numbers during the collision are the parity 
p = (—1)"+^+! and the total projection quantum num- 
ber M = rrij + Ml ■ The calculations in the present work 
are performed for p = —1 and M = 1. This choice in- 
cludes s-wave scattering of NH molecules in initial state 
nij = +1, which is magnetically trappable, to nij = 
and —1, which are not. The basis set used included all 
functions up to rimax = 1 and Lmax = 3. This uncon- 
verged basis set is sufficient for the purpose of comparing 
MQDT results with full coupled-channel calculations. 



A. Numerical methods 

The coupled-channel calculations required for both 
MQDT and the full coupled-channel ap pro ach were car- 
ried out using the MOLSCAT package [5311. as modified 
to handle collisions in magnetic fields [50] ■ The cou- 
pled equations were solved numerically using the hybrid 
log-derivative propagator of Alexander and Manolopou- 
los [53|, which uses a fixed-step-size log-derivative prop- 
agator in the short-range region (rmin < r < r^i^) and 
a variable-step-size Airy propagator in long-range region 
(''mid ^ r < r,„ax. The full coupled-channel calculations 



used Tmin — 2.5 A, r,„id = 50 A and r^ax = 250 A (where 
1 A = 10~^° m). MQDT requires coupled-channel calcu- 
lations only from rmin to Tmatch (which is less than rmid), 
so only the fixed-step-size propagator was used in this 
case. 

The MQDT reference functions and quantum defect 




FIG. 1: (Color online) Zero-field reference potentials. For the 
Vo reference potential the first rotational excited state is also 
shown (n — 1). The hard wall at r = 4.5 A is shown as a 
dashed line. 



parameters were obtained as described in Section FlI C[ us- 
ing the Numerov propagator [ssj to solve the I-dimension 
Schrodinger equations. Use of the renormalized Nu- 
merov method [56j was not found to be necessary in the 
present case. The MQDT Y matrix was then obtained 
by matching to the log-derivative matrix extracted from 
the coupled-channel propagation at a distance rmatch- 



Comparison of full coupled-channel and MQDT 
results 



1. Choice o/rniatch and reference potential 

One of the goals of MQDT is to obtain a matrix 
Y{E, B) in such a way that it is only weakly dependent 
on energy E and magnetic field B. However, the actual 
form of Y{E,B) is strongly dependent on the distance 
at which it is defined and the reference potentials used. 
In the present work we consider three different reference 
potentials, as shown in Figure [TJ First we define a ref- 
erence potential containing a pure Cg long-range term, 
which has been used with great success in cold atom- 
atom collisions. 



■ (0 = T + 



6 



E? 



(29) 



where Cg = 7.621 x 10^ A cm"! for Mg-fNH [51]. Sec- 
ondly we define a reference potential containing an addi- 
tional Cs term, 

W = + ^2 ' (30) 

where Cg = 9.941 x 10^ A cm'^ J51j. Finally we define 



[/-f^^o(r) = Vo(r) 



h'L,(L, + l) 
2/ir^ 



E° 



(31) 
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Collision Energy (K) 

FIG. 2: (Color online) The squares of diagonal T-matrix el- 
ements in the incoming channels for rrij = -f 1, L = and 2 
at -B = 10 G, obtained from full coupled-channel calculations 
(solid, black) and MQDT with the Ce reference potential and 
fmatch = 20 A (dashed, red). T-matrix elements are labeled 
with quantum numbers rrij , L, Ml- [Units of Gauss rather 
than Tesla, the accepted SI unit of magnetic field, have been 
used in this paper to conform to the conventional usage of 
this field.] 



where Vo{r) is the isotropic part of the interaction po- 
tential, which is equivalent to the diagonal W matrix 
element in the incoming s-wave channel. Each refer- 
ence potential contains a hard wall at r = rf^", so that 
U![''^{r) = cx) for r < rf'^". This allows the phase of 
the reference functions in each channel to be adjusted if 
required. A useful feature of MQDT, to be explored in 
future work, is that the position of the hard wall can be 
chosen to minimize the energy- dependence of Y . How- 
ever, in the present paper we simply take r^^" = 4.5 A. 

It is convenient to compare MQDT and coupled- 
channel results at the level of T-matrix elements, Ty = 
Sij — Sij . In general we label elements Ta^L,ML-^a',L',Mi j 
where |a) = \nsjmj). However, the collisions considered 
in the present paper are all among the n = 0, j = 1 levels 
and so a is simply abbreviated to rrij . The spin-changing 
cross sections are quite small except near resonances, so 
we focus mostly on diagonal elements, for which we sup- 
press the second set of labels. 

Figure [2] compares diagonal T-matrix elements |Tj,ip 
obtained from full coupled-channel calculations with 
those from the MQDT method for the pure Cq refer- 
ence potential of equation (I29p . with a matching dis- 
tance of Tniatch — 20 A. The Y matrix was recalculated 
at every energy at which full coupled-channel calcula- 
tions were performed. The MQDT results reproduce the 
coupled-channel results almost exactly at collision ener- 
gies E/k-Q > 10 mK. However, at lower energies the re- 
sults start to differ noticeably. It may be noted that 

prei,Ce _ „ Q Q at Wch = 20 A. 

Figure [3] shows the diagonal Y elements correspond- 
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FIG. 3: (Color online) Diagonal Y matrix elements as a func- 
tion of collision energy at B = 10 G for the Ce reference po- 
tential with r-jnatch = 20 A. The dashed vertical lines show the 
positions of quasibound states as described in the text. 



ing to Figure [2l They vary smoothly across most of the 
energy range, and are continuous across the threshold at 
zero energy, but exhibit occasional sharp structures as 
a function of energy. These sharp features are close to 
the energies of quasibound states, as shown by carrying 
out bound-state calculations using the BOUND package 
[57} . with the same basis set as the MOLSCAT calcula- 
tions. The resulting bound-state energies are shown in 
Figure |3] as dashed vertical lines. The broad feature near 
E/k^ = 0.5 K is due to a quasibound state (Feshbach 
resonance) with quantum numbers n = 1, j = 0, rrij = 0, 
i = 3. 

For MQDT to be more efficient than full coupled- 
channel calculations, it needs to produce results in agree- 
ment with full coupled-channel calculations from an 
energy-insensitive Y matrix that can be assumed to be 
constant or can be obtained by interpolation from a few 
energies, instead of being recalculated at every energy. 
However, the Y matrix elements in Figure [3] do not meet 
this requirement: the resonant features prevent reliable 
interpolation over useful ranges of energy. 

The energy sensitivity of the Y matrix in Figure [3] 
is due to the value used for r„iatch- When rmatch is 
large, resonance features due to quasibound states may 
be present in the log-derivative matrix from which Y 
is obtained. In this case the open and closed-channel 
blocks of Y are uncoupled, so that Y l^o, and the 
resonances appear through the Yqo term in Eg. 12 II rather 
than through tani' + l^c [58]. However, if Tmatch is small 
enough, the resonance features are shifted to high en- 
ergies, out of the region of interest. It is usually desir- 
able to obtain 1^ at a value of rmatch that is in or near 
the classically allowed region for all weakly closed chan- 
nels. However it must be remembered that the MQDT 
method neglects interchannel couplings that occur out- 
side Tmatch, SO there is always a tradeoff between choosing 



7 




Collision Energy (K) Collision Energy (K) 



FIG. 4: (Color online) The squares of diagonal T-matrix el- 
ements in the incoming channels for rrij = +1, L = and 2 
at -B = 10 G, obtained from full coupled-channel calculations 
(solid, black) and MQDT with the Ce reference potential and 
rmatch = 6.8 A (dashed, red). 



FIG. 5: (Color online) The squares of diagonal T-matrix el- 
ements in the incoming channels for rrij — +1, L — and 2 
at B = 10 G, obtained from full coupled-channel calculations 
(solid, black) and MQDT with the Ce + Cs (dot-dashed, blue) 
and Vo (dashed, red) reference potentials and rmatch = 6.8 A. 



a value that minimizes the energy-dependence and one 
that takes account of coupling at relatively long range. 
This is particularly important in molecular scattering, 
where the anisotropy of the interaction potential often 
provides substantial couplings at long range. 

It is convenient to consider lengths and energies in ul- 
tracold scattering in relation to the Van der Waals char- 
acteristic length and energy, defined by (soj 



''VdW 



^ 1 f 2fiCe 

2 I 



and Eydw 



(32) 



For Mg+NH, rydw = 12.7 A and -Bydw/^B = H niK. 
In atomic systems, it is common to place Tmatch close to 
Tydw- However, the quasibound state responsible for the 
broad feature in Figure [3] is due to an n = 1 state, with 
an outer turning point around 5.7 A. The resonant fea- 
ture therefore does not shift in energy significantly until 
''match is around 7 A. In addition, it is not enough simply 
to move rniatch to short range with the same reference 
function. Figure H] shows diagonal T-matrix elements ob- 
tained by MQDT with the Cg reference function, as in 
Figure [21 but with Tmatch = 6.8 A. This does indeed pro- 
duce a Y matrix without poles in the energy region of 
interest, but the MQDT results are no longer in agree- 
ment with the full coupled-channel results at any of the 
energies considered. This is because the difference be- 
tween the reference potential and the diagonal W matrix 
elements at rmatch = 6.8 A is \U-^^'^'' — Wii\/k^ « 4 K, as 
seen in Figure [TJ Alternatively, in terms of the approach 
of Mies and Raoult [13], 6.8 A is too short a distance for 
the y matrix evaluated with the Cg reference potential 
to have reached its asymptotic value. 

This problem may be remedied by using a better ref- 
erence potential. Figure [S] shows results obtained using 
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FIG. 6: Ratio of the square of the diagonal T matrix element 



for m, 



-fl, L = at B = 10 G for MQDT, with the 



Vb reference potential and rnmtch 
coupled-channel calculations. 



6.8 A, to that from full 



the reference potentials of equations (|30l) and (PT|) . again 
for rmatch = 6.8 A. The Cg + Cs reference potential gives 
a marked improvement over the pure Cg reference po- 
tential. The T matrix elements it produces follow the 
form of the full coupled-channel results but still become 
poor at energies much below 1 K: at rmatch = 6.8 A, 
jj^rof,C6+C8 _ 0.35 K. However, the resuhs 

obtained with the Vq reference potential are much more 
accurate, and can scarcely be distinguished from the full 
coupled-channel results in Figure [5] 

Even the Vb reference potential does not produce ex- 
act results. Figure [S] shows the ratio of the MQDT T- 
matrix elements for this reference potential to the full 
coupled-channel results. The poles in the ratio arise sim- 
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FIG. 7: (Color online) Diagonal Y matrix elements as a func- 
tion of energy at B = 10 G, for the Vo reference potential with 
r-match = 6.8 A. 
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FIG. 8: (Color online) Diagonal Y matrix elements as a func- 
tion of magnetic field at E/ks = 400 mK, for the Vo reference 
potential with rmatch = 6.8 A. 



ply because MQDT places the zeroes in |T|2 (where the 
phase shift is an integer multiple of tt) at very slightly dif- 
ferent collision energies. However, at very low energies 
(below about 1 mK) the MQDT results underestimate 
the squared T-matrix elements by up to 3%. This prob- 
ably arises because the "best" reference potential would 
be one that takes account of adiabatic shifts due to mix- 
ing in excited rotational levels. For the n = channels, 
the shift due to n = 1 channels may be estimated from 
2nd-order perturbation theory to be about 0.012 cm~^ 
(equivalent to 17 mK) at rmatch = 6.8 A. This will cause 
residual errors in the MQDT C functions that are re- 
sponsible for the small errors visible in Figure [B] 

Figure [7] shows representative matrix elements of Y 
obtained at rmatch = 6.8 A, with the Vq reference poten- 
tial, as a function of energy. It may be seen that they are 
nearly linear in energy. The other matrix elements of Y 
show similar behavior. While the actual values of matrix 
elements vary substantially, they are all nearly linear in 
energy for rmatch = 6.8 A. 

It should be noted that when the reference functions 
are obtained numerically, as in the present work, there 
is no significant difference in computer time for different 
choices of reference potential. Using the full Vq reference 
potential is just as inexpensive as using a simpler one. 



2. Feshbach resonances 

Magnetic fields have important effects on cold molec- 
ular collisions, and in particular magnetically tunable 
low-energy Feshbach resonances provide mechanisms by 
which the collisions may be controlled. It is therefore 
important to establish whether the Y matrices obtained 
from MQDT are smooth functions of magnetic field as 
well as energy and can be used to characterize Feshbach 
resonances. If they are, it will offer substantial compu- 



tational efficiencies. 

Figure H] shows how the diagonal Y matrix elements 
vary as a function of magnetic field for Mg-|-NH collisions 
over the range from to 2500 G for a collision energy of 
400 mK. It may be seen that the matrix elements are 
indeed very nearly linear, as required for efficient inter- 
polation. 

In Mg-fNH, there is a Feshbach resonance due to the 
n ~ I, j — 0, rrij = 0, i = 3 state shown in Figure [3] 
that tunes down towards the n — 0, ruj = -1-1 threshold 
with increasing field. Figure ^ shows the comparison be- 
tween MQDT and full coupled-channel calculations for a 
selection of diagonal and off-diagonal T matrix elements 
as the magnetic field is tuned across this resonance at 
energies of 400 mK and 1 mK. At each energy, MQDT 
results were obtained both by recalculating the Y ma- 
trix at every field and by linear interpolation between 
two points separated by 100 G. In both cases, the inter- 
polated MQDT results are within about 0.2 G of the full 
MQDT results even for this long interpolation, and this 
could of course be improved simply by considering a few 
more fields across the range to allow. However, there is 
also a residual error of of 0.1 to 0.2 G in the resonance po- 
sition even for the full MQDT results, which is not very 
different at the two collision energies considered. This is 
again likely to be due to the effect as described in Sec- 
tion [TTrBTJ the Vq reference potential neglects couplings 
between channels outside rmatch, and for the small value 
of rmatch uscd here these couplings are sufficient to shift 
the resonance positions slightly. Apart from these small 
shifts, however, both the elastic and the inelastic scat- 
tering around the resonances are very well described at 
both energies. 

The linearity of the Y matrix with both energy and 
applied magnetic field is an extremely promising result, 
and suggests that MQDT will provide very efficient ways 
of performing cold collision calculations as a function of 
energy and magnetic field, without needing to repeat the 
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FIG. 9: (Color online) Squares of diagonal and off-diagonal 
T-matrix elements as the field is tuned across a Feshbach 
resonance at E/k^ = 400 mK (upper panel) and 1 mK (lower 
panel). The MQDT results are obtained with the Vb reference 
potential at rmatch = 6.8 A. 

expensive coupled-channel part of the calculation on a 
fine grid. 

IV. CONCLUSIONS 

We have shown that Multichannel Quantum Defect 
Theory (MQDT) can be applied to low-energy molecular 
collisions in applied magnetic fields. MQDT provides a 
matrix Y , defined at a distance rmatch at relatively short 
range, which encapsulates all the short-range dynamics of 
the system. For the prototype Mg-|-NH system, we have 
shown that MQDT can provide numerical results that 



are in quantitative agreement with full coupled-channel 
calculations if the MQDT reference functions are defined 
appropriately. 

We have investigated the effect of different choices of 
reference potential and values of rmatch- For cold atom- 
molecule collisions, unlike cold atom-atom collisions, cal- 
culations are likely to be needed over a significant range 
of collision energy, perhaps 1 K or so. If rmatch is placed 
at too long a range, there is a significant likelihood of 
resonant features within the energy range that prevent 
simple interpolation of Y. This may be circumvented by 
carrying out the matching at a smaller distance rmatch- 
However, when this is done, a pure Cq reference poten- 
tial may not be sufficient. For Mg-|-NH, the most satis- 
factory procedure is to perform matching at fairly short 
range (inside 7 A) and use a reference potential that is 
defined to be the same as the true diagonal potential in 
the incoming channel. 

The major strength of MQDT for molecular applica- 
tions is that, if the the matching to obtain Y is carried 
out at relatively short range, the matrix is only weakly 
dependent on collision energy and magnetic field. This 
allows very considerable computational efficiencies, be- 
cause the expensive calculation to obtain Y needs to be 
carried out at only one or a few combinations of colli- 
sion energy and field. The remaining calculations to ob- 
tain scattering properties on a fine grid of energies and 
fields are then computationally inexpensive, varying only 
linearly with the number of channels N. Full coupled- 
channel calculations, by contrast, scale as N^. 

MQDT is a promising alternative to full coupled- 
channel calculations for cold atom-molecule collisions, 
particularly when fine scans over collision energy and 
magnetic field are required. In future work, we will in- 
vestigate further the choice of reference functions to op- 
timize the accuracy and to minimize the dependence of 
Y on collision energy and field. We will also investigate 
how the results for Mg-I-NH transfer to more strongly 
anisotropic systems, with stronger long-range anisotropy 
and more closed channels that are capable of producing 
scattering resonances. 
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